Video-rate laser Doppler vibrometry by heterodyne holography 
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We report a demonstration video-rate heterodyne holography in off-axis configuration. Reconstruction and 
display of 1 Megapixel holograms is achieved at 24 frames per second, with a graphics processing unit. Our 
claims are validated with real-time screening of steady-state vibration amplitudes in a wide-field, non-contact 
vibrometry experiment. © 2013 Optical Society of America 
OCIS codes: 090.1995, 280.3340 
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The laser Doppler method is the most common optical 
interferometry technique used for non-contact measure- 
ments of mechanical vibrations. Though highly effec- 
tive for single-point vibration analysis, this technique 
is much less adapted to wide-field imaging than holog- 
raphy. Homodyne [IH3] and heterodyne [3H6] holo- 
graphic recordings in off-axis configuration enabled re- 
liable measurements of mechanical vibrations, but none 
of them allowed real-time monitoring, which is an es- 
sential feature. Matching the display rate of optically- 
measured Megapixel digital holograms with real-time 
imaging standards is demanding in terms of computa- 
tional power. Holographic measurements are performed 
in a diffraction plane. Hence, image formation requires to 
simulate the back-propagation of an optical field. Such 
propagation involves turning the data measured in the 
plane of an array detector into a reciprocal plane with at 
least one bi-dimensional numerical Fourier transforma- 
tion, typically a Fast-Fourier transform (FFT). Recently, 
real-time display of digital holograms with Graphics Pro- 
cessing Units (GPUs) [TIE] has alleviated the issue of 
the high computational workload needed for such im- 
age reconstruction. Parallel computations on GPU con- 
sistently increase the throughput with respect to CPU 
for computer-generated holograms, which demonstrated 
the performance of GPUs in streamline image process- 
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In this letter, we report an experimental demonstra- 
tion of an image acquisition scheme designed to perform 
video-rate image reconstruction and display from hetero- 
dyne holographic measurements on a 1 Megapixel sen- 
sor array. Image reconstruction of steady-state vibration 
modes up to 100 kHz at a rate of 24 images per second is 
achieved. GPU processing is shown to enable holographic 
reconstruction and display with 3 FFT calculations per 
recorded frame, which covers the processing throughput 
needs of three reconstruction approaches : the convo- 
lution, angular spectrum, and Fresnel transform meth- 
ods p]. 

The acquisition setup consists of an off-axis, 
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Fig. 1. Experimental image acquisition setup. 



frequency-shifting holographic scheme, used to perform 
a multipixel heterodyne detection of optical modulation 
sidebands. Optical heterodyning is a process for plac- 
ing information at frequencies of interest (e.g. the me- 
chanical vibration of an object under investigation) into 
a useful frequency range by mixing the frequency con- 
tent of the probe beam with a reference (or local os- 
cillator, LO) beam. The optical frequency of the refer- 
ence beam is shifted to generate a beat frequency of 
the interference pattern within the sensor bandwidth, 
which carries the information at the original frequency 
of interest. The Mach-Zehnder heterodyne interferom- 
eter used for the detection of an object field E in re- 
flective geometry, beating against a LO field E^o, is 
sketched in fig. [TJ The main optical radiation field is pro- 
vided by a 100 mW, single-mode, doubled NdiYAG laser 
(Oxxius SLIM 532) at wavelength A = 532 nm, and op- 
tical frequency cjl/(27t) = 5.6 x 10 14 Hz. The optical fre- 
quency of the LO beam is shifted by an arbitrary quan- 
tity Auj in the low radiofrequency (RF) range (0 Hz < 
Acj/(27r) < 100 kHz) by two acousto-optic modulators 
(AA-electronics, MT80-A1.5-VIS). The LO field has the 
form £?lo = ^log^ u;l+Au; ^, where £lo is its complex 
amplitude. The object studied is a sheet of paper, whose 
lateral dimensions are 9 x 26 mm, shined over 9x17 mm. 
It is attached to a piezo-electric transducer (PZT, Thor- 
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labs AE0505D08F), vibrating sinusoidally, driven at 10 
V. Its local vibration is z(t) = z max sin(o;M^), where z max 
is the vibration amplitude, provokes a modulation of the 
optical path length of the object field. It induces a local 
optical phase modulation of the backscattered field at 
the mechanical vibration frequency <j)(t) = </>o sui(co>m£), 
where 0o = 47rz max /A is the modulation depth of the 
optical phase. Holographic analysis of objects in sinu- 
soidal vibration with a frequency-shifted LO beam tuned 
to the n-th optical modulation sideband was carried-out 
extensively in [4;. The frequency filtering properties of 
time-averaged holography were introduced in pQ. The 
same filtering properties were described in digital holog- 
raphy [l2l[T3]. It was shown that holograms of optical 
sidebands could be recorded selectively by matching the 
frequency shift of the local oscillator beam with the fre- 
quency of the band of interest. In the case of sinusoidal 
phase modulation, the hologram amplitude is of the form 
Jn(4>o), where J n is n-th order Bessel function of first 
kind. This modulation yields fringes that correspond to 
local extrema of J n . It was also shown that imaging at 
harmonics of the vibration frequency could enable ro- 
bust assessment of vibration amplitudes which are much 
greater than the optical wavelength [6] . The interference 
pattern / is measured with a Pike F421-B camera on a 
Kodak KAI-04022 interline, progressive-scan CCD sen- 
sor (2048 x 2048 pixels, pixel size d px = 7.4 /im). The 
camera is run in binning mode (an effective pixel is made 
of 4 adjacent pixels); 16 bit, 1024 x 1024 pixels images 
are sampled at o;s/(27r) = 24 Hz throughout the experi- 
ments described hereafter. The RF command signals at 
frequencies Acj, cjm, and ujs are phase- locked. The tem- 
poral part of object field undergoing sinusoidal phase 
modulation can be decomposed in a basis of Bessel func- 
tions using the Jacobi- Anger identity 



E = g e «^t+*0(t) = J2 SJ n (Me 
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where £J n ((j)o) = S n is the weight of the optical mod- 
ulation sideband of order n, and where £ is the com- 
plex amplitude of the field. If the frequency detuning 
Auj is set close to the n-th modulation harmonic, i.e. 
\Auj — uujm\ < <^s? and if the modulation frequency is 
much greater than the sampling frequency, i.e. ujm ^> ^s? 
the time-averaging-induced bandpass filter of the detec- 
tion process will isolate the term of order n in eq. [T] and 
reject all other optical sidebands. In the sensor plane, 
the interference pattern of E and £?lo takes the form 
I(t) = \E+ E LO \ 2 = \E\ 2 + \E LO \ 2 + EEl Q + E*E LO , 
where * denotes the complex conjugate. The frame I(t) 
acquired at time t by the framegrabber is moved to a 
frame buffer in the GPU RAM by a CPU thread (fig.EJ). 

To detect the heterodyne signal of interest EE^ Q = 
S^ S n ^ nuJm ~ - Auj ) f \ a sliding 4-phase temporal demodu- 
lation is performed. The intermediate frequency uujm — 
Auj is set within the camera bandwidth to be sam- 
pled efficiently. More precisely, for a LO detuning Auj = 
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Fig. 2. Algorithmic layout of holographic rendering. 

tiujm — cjs/4, the modulation sideband £ n beats at the 
frequency cjs/4 (6 Hz) in I(t). To detect it, the following 
quantity is formed 



H(t) = ^2l(t-2 P 7r/uJs)e^^ 2 



(2) 



p=0 



H is a complex-valued array, proportional to the 
diffracted field E in the sensor plane. Its calculation re- 
quires the allocation of a stack of 4 arrays in the GPU 
memory, filled with /(£), I(t — 27t/cjs), I(t — 47t/cjs), 
and I(t — 6tt/ujs). Each new frame grabbed at instant t 
yields a shift of the stack (fig.EJ) : I(t) replaces the array 
J(t-87r/a;s). 

In off-axis configuration, the spatial spectrum in the 
reciprocal plane (k x ,k y ) of the term EE^ Q in the ex- 
pression of I(t) is shifted by the projection of the wave 
vector difference along a transverse direction x in the 
sensor plane Ak x ~ 2i7r#o/A. Spatial filtering [14] of 
the time-demodulated signal H(t) in off- axis geometry is 
used to remove the remaining contributions of the zero- 
order terms \E\ 2 and |£^lo| 2 and the twin-image term 
E*E\,o to enhance the detection sensitivity. It is made 
by multiplying the (k Xl k y ) spectrum of H by a mask M, 
allowing only frequencies in the neighborhood of Ak to 
pass. 

F(t) = F- 1 {MF{H{t)}} (3) 

where J 7 is a spatial FFT and T~ x is an inverse spa- 
tial FFT. This operation is handled by the GPU. Only 
the heterodyne contribution of interest F(t) = £l £ n 
remains in the filtered frame. 

Image rendering from F oc S n in the sensor plane in- 
volves a scalar diffraction calculation in the Fresnel ap- 
proximation, performed with a discrete Fresnel trans- 
form [15] . The hologram f(t) back-propagated to the 
object plane is calculated by forming the FFT of the 
product of F with a quadratic phase map, depending on 
the relative curvature of the wavefronts of E and £?lo 
in the sensor plane (x, y) via a distance parameter Az. 
This calculation is handled by the GPU. 



f(t)=?{F{t)e 



} 



(4) 



Finally, the GPU calculates the quantity \f(t)\ 2 oc 
| J n (47rz max /A)| 2 , which is a map in the object plane of 
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(b) 
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the composition of the local vibration amplitude field 
z max with the squared amplitude of the Bessel function of 
order n. Image brightness adjustment is also performed 
by the GPU. Those maps are displayed in fi^3] for the 
modulation sidebands of order n = (a), n = 1 (b), 
n = 3 (c), and n = 7 (d). Excited at cjm/(27t) = 10 kHz, 
the paper sheet builds up a steady-state vibrational 
mode with rectilinear nodes and bellies oriented along y, 
with a mechanical wavelength of ~ 5 mm. In fig.[3ja), the 
non-moving support of the object in vibration is visible 
(arrow). A sweep of the detection sideband is reported in 
media 1. Additional vibrational patterns are screened for 
excitation frequencies swept from Hz to 20 kHz (media 
2) and from Hz to 100 kHz (media 3), with a detection 
tuned to the first modulation sideband. The propensity 
of the filter of eq. [3] to cancel-out spurious artefacts is 
assessed in real time in media 4. 

The image reconstruction and display algorithm was 
elaborated with Microsoft Visual C++ 2008 integrated 
development environment and NVIDIA's Compute Uni- 
fied Device Architecture (CUD A) software development 
kit 3.2. FFT calculations were made with the function 
cufftExecC2C() from the CUFFT 3.2 library on single 
precision floating point arrays. The program was com- 
piled and run on Microsoft Windows 7-64 bit. The 
computer hardware configuration was based on an ASUS 
P6T motherboard with a 2.67 GHz Intel core i7 920 CPU 
and a NVIDIA GeForce GTX 470 GPU. Image rendering 
calculations I ^ H ^ F ^ \ f\ 2 are performed sequen- 
tially, in the main GPU thread (fig. [2]); the whole process- 
ing time of one frame is reported in table [TJ For bench- 
mark purposes, the rendering performance of 2048 x 2048 
pixels recordings read out at cos / (2tt) = 8 Hz is also re- 
ported in table [TJ 



Table 1. Benchmarks of image rendering time. 



Array size (pixels) 


time 


1024 x 1024 


5.5 - 5.9 ms 


2048 x 2048 


19.7 - 20.4 ms 



We have demonstrated that the detection and ren- 
dering of 1 Mega Pixel heterodyne holograms can be 
carried-out with a refreshment rate of 24 Hz with com- 
modity computer graphics hardware. Video-rate optical 
monitoring of steady-state out-of-plane vibration ampli- 
tudes was reported. This demonstration opens the way 
to high throughput multipixel optical heterodyne sens- 
ing in real-time. 

This work was funded by Fondation Pierre- Gilles 
de Gennes (FPGG014 grant), Agence Nationale de la 
Recherche (ANR-09-JCJC-0113 grant), And Region Re- 
de-France (C'Nano grant). 
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